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Abstract 

In this work we study the stability of the Jordan-Brans-Dicke (JBD) static universe. This is 
motivated by the possibility that the universe might have started out in an asymptotically JBD 
static state, in the context of the so called emergent universe scenario. We extent our previous 
results on stability of JBD static universe by considering spatially homogeneous Bianchi type IX 
anisotropic perturbation modes and by including more general perfect fluids. Contrary to general 
relativity, we have found that the JBD static universe, dominated by a standard perfect fluid, could 
be stable against isotropic and anisotropic perturbations. The implications of these results for the 
initial state of the universe and its pre-inflationary evolution are discussed. 
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I. INTRODUCTION 



Measures coming from three different independent sources (CMB, Type la supernova, and 
duster abundances) strongly suggest that the expansion of the universe has been accelerating 
during recent epoch [1]. This discovery has become one of the main challenges for modern 
theoretical physics. 

In general, dark energy or quintessence, responsible for the cosmic acceleration, deter- 
mines the features of the evolution of the universe. Certainly, the nature of this sort of 
energy may lead to the improvement of our picture about particle physics and/or gravita- 
tion. Most of these studies have been done in the standard theory of gravity, i.e. general 
relativity theory (GR). However, motivated mainly from string theories, a less standard the- 
ory have been carried out, namely the so called scalar-tensor theory ofgravityflQQQ. An 
important advantage of these models is that they naturally allow |3j, |5| a super-accelerating 



expansion of the universe where the effective dark energy equation of state w = - crosses 



the phantom divide line w = —1. Such a crossing is consistent with current cosmo logical 

data 61 and if it is confirmed, it would become an enigmatic problem which can not been 

□ fl 

explained easily with standard quintessence models [7, &] . 

The archetypical theory associated with scalar tensor models is the JBD gravity. The 
JBD theory [2j is a class of models in which the effective gravitational coupling evolves with 
time. The strength of this coupling is determined by a scalar field, the so-called JBD field, 
which tends to the value the inverse of the Newton's constant. The origin of JBD theory 
is in Mach's principle according to which the property of inertia of material bodies arises 
from their interactions with the matter distributed in the universe. In modern context, JBD 
theory appears naturally in supergravity models, Kaluza-Klein theories and in all known 



effective string actions 
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The study of static universe and its stability has always been of great interest since the 
pioneer work of Eddington [if]]. For example, in the context of GR the stability of the 



Einstein static ( 
in Refs. 



jS) universe in the presence of conventional matter field has been studied 



191 ]. In the presence of ghost scalar field it was studied in Ref. 20J. Also, 



the ES universe has been studyied in different gravitational theories. In GR a generalization 



which include a variable pressure have been ana. 
world models it was considered in Refs. 22|, |23 . 



yzed in Ref. 



24 



211 ] . In the context of brane 



25) . The study in the Einstein-Cartan 



2 



theory it is 



in Refs. 



27 



bund in Ref. 



28 



261 ] . In loop quantum cosmology, this subject has been studied 



29(. The stability of the ES universe in a f(R) gravity and in modified 



Gauss-Bonnet gravity theories have been studied in Refs. [30( and [3l|, respectively. 

Recently, the stability of ES models has become relevant for the study of cosmological sce- 



narios in which the ES universe corresponds to an initia 
cosmology, the so-called emergent universe scenario 



universe 



32 



state for a past-eternal inflationary 
. The original idea of an emergent 
33] is that in which the universe emerges from an ES universe state with radius 
a » l p (where a$ is the scale factor at some instant and l p is the Planck length), inflates 
and then is subsumed into a hot Big Bang era. Such models are appealing since they provide 
specific examples of nonsingular (geodesically complete) inflationary universes. Also, these 
models could avoid an initial quantum-gravity stage if the static radius is larger than the 
Planck length. 



However, the emergent universe models based on GR, with ordinary matter, suffer 



a number of important shortcomings. In particular, the instability of the ES state 23 



rom 



28 



29] makes it extremely difficult to maintain such a state for an infinitely long time. The 



instability of the ES solution ensures that any perturbations, no matter how small, rapidly 
force the universe away from the static state, thereby aborting the scenario. 

Some models have been proposed to solve the stability problem of the asymptotic static 
solution. They consider non-perturbative quantum corrections of the Einstein field equa- 
tions, either coming from a semiclassical state in the framework of loop quantum gravity 
(LQG) 28|, |29j or braneworld cosmology with a timelike extra dimension 23j, |24[. Other 
possibilities to consider are the Starobinsky model or exotic matter 34. l35|. 

On the other hand, it has been shown that a scalar tensor theory could solve the problem 



361 ] . it was found 



of the instability of the emergent universe models. In particular, in Ref. 
that a self interacting JBD theory presents a stable past eternal static solution, which even- 
tually enters a phase where the stability of this solution is broken leading to an inflationary 
period, providing in this way, an explicit construction of an emergent universe scenario. 

In this work, we study the stability of the JBD static universe. This is motivated by 
the possibility that the universe might have started out in an asymptotically JBD static 
state [36] . We extent our previous results on the stability of JBD static universe by consider 



spatially homogeneous Bianchi type IX anisotropic perturbation modes and by including 
more general perfect fluids. General anisotropic perturbations are important to be consider 



3 



because they could be the crucial destabilizing modes of a static universe, see Refs. Q]- 
Contrary to the GR case we have found that the JBD static universe dominated by a perfect 
fluid could be stable against isotropic and anisotropic perturbations for some sort of matter 
components. 

The paper is organized as follows. In Sect. [Ill we review briefly the cosmological equations 
of the JBD model. In Sect. HHI the existence and nature of static solutions are discussed for 
universes dominated by a standard perfect fluid and by a scalar field. In Sect. [IV] we study 
the stability of the JBD static universe against small anisotropic perturbations. In Sect. [V] 
we focus in a particular example of a JBD potential which allow us to use a dynamical 
system approach to study the problem of stability of the JBD static universe. In Sect. [VT] 
we summarize our results. 



II. THE MODEL 



We consider the following JBD action for a self-interacting potential and matter, given 



by 



S 



d x 



I $jR _ I| V ^V^ + V(<S>) + C r 



(1) 



where C m denote the Lagrangian density of the matter, R is the Ricci scalar curvature, $ is 
the JBD scalar field, w is the JBD parameter and V(<&) = V is the potential associated to 
the field $. In this theory 1/$ plays the role of the gravitational constant, which changes 
with time. This action also matches the low energy string action for w = —1 
From the Lagrangian density, Eq. ([!]), we obtain the field equations: 



15|- 



R 



w 



1 r, ^ x. 



$ 2$ 



— T 



(2) 



and 



□$ = — - — T\ + — - — 

3 + 2w M 3 + 2w 



dV 



(3) 
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where we have consider that $ is a function of the cosmological time,t, only. Units are such 
that c = h = 1. 



III. THE STATIC UNIVERSE SOLUTION IN JBD THEORY 

Let us start by considering the closed Friedmann- Robert son- Walker metric: 

dr 2 



ds 2 = dt 2 - a(ty 



+ r 2 (d6 2 + sin i 



(4) 



1 — r 2 

where a(t) is the scale factor, t represents the cosmic time. The matter content of the universe 
is modelled by a perfect fluid with effective equation of state given by P = (7 — 1) p. In 
general, when the perfect fluid is described by a scalar field, it is found that the parameter 
7 becomes variable. Thus, by using the metric, Eq. (jl]), the set of field equations (j2J) and 
(J3J) become 

TT 2 1 P w f ®Y V 



a 



and 



a rr2 1 $ T T $ w / $ \ V P 



$ + 3H$ = ) P 3P ] + [2 V - $ V] . (7) 

(2w + 3) 2w + 3 L J V 7 

The energy-momentum conservation implies that 

p + 3H(p + P) = , (8) 

where V = <iU($)/c?$. Dots mean derivatives with respect to the cosmological time. 

In the context of JBD theory the static solutions are closed universes characterized by 
ihe conditions a = a = Const., do = = do and $ = $ = Cte., $0 = = $0, see Ref. 



36]. 



Then the static solution for a universe dominated by a general perfect fluid is obtained 
if the following conditions are fulfilled 



V ' 



Po 



v 



(9) 



(10) 
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2 (. V \ n $ 



and 

7o = 1(1 + ^)2^, (II) 
3 V PoJ %Po 

where V = V{%) and V ' = (dV($)/d$) <i>=$ - These equations connect the equilibrium 
values of the scale factor and the JBD field with the energy density and the JBD potential 
at the equilibrium point. 

Note that in order to obtain a static solution we need to have a non-zero JBD potential 
with a non- vanishing derivative at the static point $ = $ - The original Brans-Dicke model 
corresponds to V(<&) = 0. However, non-zero V(<&) is better motivated and appears in many 
particle physics models. In particular, V(<&) can be chosen in such a way that $ is forced 
to settle down to a non-zero expectation value, $ — > m^/Rn, where m p = 10 19 GeV is the 
present value of the Planck mass. On the other hand, if V(<&) fixes the field $ to a non-zero 
value, then time-delay experiments place no constraints on the Brans-Dicke parameter w 



3 71 ] . In particular, if we choose the JBD potential in such a way that <£> will be stabilized 



at a constant value, let say $/, at the end of the inflationary period (see Ref. 36j] as an 



example), we can recover GR by setting <&/ = ml/87i, together with an appropriated value 



for the parameter w which will be in agreement with the solar system bound |37l . 



38] 



A. JBD static universe dominated by a standard perfect fluid 

The static solution is characterized by the Eqs. flWTTl). from which we obtain 7 > | as a 
condition for a static solution, if the JBD potential, V"($), is positive 36|. Notice that this 
means that it is not possible to have a static solution if the universe is dominated by the 
cosmological constant (corresponding to 7 = in the equation of state), but it is possible to 
have a static universe when it is dominated by dust or radiation, among others possibilities. 

Now, we study the stability of this solution against small homogeneous and isotropic 
perturbations. In order to do this, we consider small perturbations around the static solution 
for the scale factor and the JBD field. We set 

a(t) = ao[l + £(*)], (12) 

and 

$(t) = $ [l +£(*))]. (13) 



Then, we have 

p = p + Sp(e) w p - 37p £, (14) 

where e C 1 and (3 <C 1 are small perturbations. By introducing expressions ffl~2l) . f|T3|) and 
(fl4"|) into Eq. (jSJ) and Eq. ([7]), and retaining terms at the linear order on e and (3, we obtain 



the following coupled equations 



e — 



a, 







(f 







--£-1 = 

2 u , 




2 a2 



and 



(3 + 2w)^- ( -^-2$ V ")/3+(4-3 7 )-|^ = 



a / a 



(15) 



(16) 



where Vtf' = (d?V($)/d$% = * . 

From the system of Eqs. (ITS"]) and ffTB"]) we can obtain the frequencies for small oscillations 



a 2 Q (3 + 2w) 



a2$ y "-6 + w(2-3 7 ) 



(17) 



± yj [ - 6 + al<& V£ + 2w - 3w 7 ] +2(3 + 2w) ( - 6 + a2$ y o " [3 7 - 2]) 

Note that the static solution is stable if the inequality, u 2 . > 0, is fulfilled. Assuming that the 
parameter w satisfies the constraint, (3 + 2u>) > 0, it is found that the following inequalities 
must be achieved in order to have a stable static solution 



- < 7 < - 

3 ' 3 



or - < 7 , 
3 



(18) 



- < w < 
2 



(7-1) 
(2 - 3 7 )2 



and 



2(6 + w) -3(3 + w) 7 + v / 3|4-3 7 |v / 3 + 2w < a 2 $ Vq < 



6 



(19) 



(20) 



3 7 -2 

From these inequalities we can conclude that for a universe dominated by a standard perfect 
fluid (with 7 > 2/3), it is possible to find a solution where the universe is static and stable. 

Here, the only exception is radiation, where 7 = 4/3, which becomes explicitly excluded 
by the latter inequalities. This peculiar behavior for radiation could be understood due 
to the particular way in which the perfect fluid appears in the equation for the JBD field, 



Eq. (jTj), where it becomes independent of the energy density and pressure of the fluid. Then, 
we can note that in the radiation case w± are both real numbers. On the other hand, w\ 
could be a positive number, but w 2 ^ is always negative. Therefore, in this case, we have a 
saddle instability. 



B. Scalar fields 



In the context of emergent universe models, the static JBD universe dominated by a scalar 



361 ] . Here, we reproduce the main results concerning this 



field (inflaton) was studied in Ref. 
static solution. 

The energy density, p, and the pressure, P, are expressed by the following equations 

p =^.+u(n (21) 

and 

P = Y~ uw ■ (22) 

Here, U(^>) represents the scalar potential associated to the scalar field 

We could write an effective equation of state for the scalar field, expressed by the 
equation P — (7 — 1) p, where the equation of state "parameter", 7, could be written as 

7 = 2fl-^V (23) 



p J 

During the static regimen, in the context of an emergent universe models, the matter 
potential U(^f) is consider as a flat potential, that is U(ty) = Uq = Const, and the scalar 
field rolls along this potential with a constant velocity \EV The conditions for static universes, 
Eqs. flUl TTTj) . imply that the following condition for the state parameter 

7„ = 2-^-2 fl-^V (24, 



«o Po V Po 
must be satisfied. 

The velocity when the scalar field ^ is rolling along a constant potential, Uq, it becomes 
expressed in terms of the static values of the scale factor, a , and the JBD field, $ - It 
results to be 

^ = 2$- (25) 



s 



Assuming (3 + 2w) > 0, the following stability conditions were obtained 



0<ag$ o V r o "<|, (26) 

and 



■- < w < — - 

2 4 



'9 - 6ag$ V " + (3 + al^') . (27) 



In relation to the conditions fl26|) and (1271) let us mention that the first inequality imposes 
a condition on the JBD potential, specifically for its first and second derivatives: < Vq < 
Vq/ (2$o)- The second inequality restricts the values of the JBD parameter. Notice that this 
inequality imposes that w < 0. JBD models with negative values of w have been considered 



in the context of late acceleration expansion of the universe [391,1401, but also appear in low 
energy limits of string theory |15j]. On the other hand, as was mentioned above, we choose 
the JBD potential, V($), in such a way that $ will be stabilized at a constant value, namely 
$/ = ml/8-K. 

Thus, from Eqs. fl26|) and (T27j) we can conclude that for a universe dominated by a 
scalar field it is possible to obtain a static solution, stable under homogenous and isotropic 
perturbation. 

IV. ANISOTROPIC PERTURBATIONS 

If we are interested in studying the stability of the static universe an important point, 



showed in Ref. 19|, is that the crucial destabilizing modes are not only the conformal per- 
turbations considered in the previous section. Anisotropic perturbations could be even more 
important. For example, in the case of the ES universe it is known that the static solution is 
neutrally stable to inhomogeneous scalar perturbations with high enough sound speed and 
to vector and tensor isotropic perturbations 3- However, this analysis does not cover 
spatially homogeneous, but anisotropic modes. It turns out that there are various unstable 
spatially homogeneous anisotropic modes [19|. This suggest that anisotropic perturbations 
could be the crucial destabilizing modes. 

In this section we proceed to study the stability of the static solution found in the previous 
section against theses anisotropic perturbations modes. In particular, we consider the general 
case of spatially homogeneous Bianchi type IX perturbations modes. 



9 



In this context, the JBD static universe is a particular exact solution of the Bianchi type 
IX, or Mixmaster universe, containing a perfect fluid. The Mixmaster is a spatially homo- 
geneous closed (compact space sections) universe of the most general type. It contains the 
closed isotropic Friedmann universes as particular cases when a fluid is present. Physically, 
the Mixmaster universe arises from the addition of expansion anisotropy and 3-curvature 
anisotropy to the Friedmann universe. 

The diagonal type IX universe has three expansion scale factors, i.e., a(t), b(t) and c(t), 
and the diagonal Bianchi IX metric is expressed by 



ds 2 



dt 2 - r] af3 (t)w a w p 



(28) 



where 



\ 



a 2 (t) 

V*f3(t)= b 2 (t) , (29) 
\ c 2 (t) / 
and the w a are differential 1-forms invariant under 5*0(3) transformation. 

In the following we will consider a universe dominated by a general perfect fluid whose 
equation of state is P = (7 — l)p. We will assume that $ and p are function of the time t, 
only. Then, by using the metric, Eq. (|28|) . in the action ([1]), we obtain the following set of 
equations for the non-null components. The (0,0) component becomes 



a 2 b 2 c 2 ab at be 

la 2 ' lb 2 ' 2c 2 46 2 c 2 4c 2 a 2 4a 2 6 2 ab ac be 



1 1 1 

+ ^ + 



P_ 



a b c \ $ w 
a b c / $ 2 



The (1,1) component is given by 




(30) 



1 1 1 3 a 2 lb 2 1 c 2 be b c 
2a 2 lb 2 2c 2 4 6 2 c 2 4a 2 c 2 4a 2 6 2 be b c 




+ 



(31) 
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The other two nonzero equations, components (2,2) and (3,3), are just cyclic changes in the 
scale factors (a,b,c) in Eq. fl3Tl) . 

The equation for the JBD field, Eq. ([3]), becomes given by 



a 
a 



b c 

1 + 

c 



* - 3P ) 



3 + 2w 3 + 2w 



2V-Q 



dV 



(32) 



On the other hand, the conservation of energy-momentum implies that 



. a b c . 

, I- + - + -K„ + p) 



0. 



(33) 



The static solution discussed in the previous section correspond to the case where a{t) = 
b(t) = c{t) = | a and $ = $ - Here, the constant values, a and $o, satisfy the conditions 
for a static solution. This was discussed in Sect. [TTJ 

In order to study the stability of this solution against anisotropic Bianchi type IX per- 
turbations, we take small perturbations around the static solutions of the scale factors and 
the JBD field. We set 



a(t) = ^[l +ei (t)l (34) 

bit) = | [l+e 2 (t)}, (35) 

and 

c(t) = j[l + es(t)], (36) 

together with the perturbation associated to the JBD field, expressed by Eq. (Tl3l) . Here, the 
parameters Ei{i = 1,2,3), just like the parameter (3, are small perturbations. Therefore, 
they satisfy <ti 1. 

For the energy density, pressure and state parameter we take 



p = p + 8p(e 1 ,£ 2 ,e 3 ) 
P = P + 6P(e 1 ,e 2 ,e 3 ) 



(37) 
(38) 



and 



7 = 7o + h( £ i, £ 2,£3) 



(39) 
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respectively. The specific form of dp, 5P and £7 depend on the kind of perfect fluid under 
consideration. 

Now, introducing these latter expressions into Eqs. (I30H33P and retaining only the linear 
terms in the perturbation parameters, we obtain that 



Q^ + 2^ + 2^ + e 2 + e 3 = -^- + 2^-(3, (40) 

- 6 £§ + 2% + 2% + 6 3 + 6 1 = ~ + 2^-'$, (41) 
a a a <p Q a 

-6% + 2% + 2% + 6\ +6 2 - S -^ + 2^-P, (42) 



and 



(3 + 2w) Q '$ = (4 - 3 7o ) Sp - 3(5 7 po + 2 0o /? - 2 <g V " (3 . (43) 
In the next subsections we study universes dominated by different type of perfect fluids. 



A. Standard perfect fluid 

Here, we consider the case of a universe dominated by a standard perfect fluid, where 7 
is a constant. Then, Eqs. ( l37|) and ( l38l) become 



and 



p = p + 5p(ei, 6 2 , 63) ~ Po - 7 Po [61 + £2 + 63] 



P = P + 5P{e u 6 2 , 63) « P + 7 (1 - 7) Po [61 + 6 2 + 63] 



(44) 



(45) 



respectively. By introducing these expressions into Eqs. ( 14011431) . and retaining the linear 
order in the parameters e»(i = 1,2, 3) and (3, we obtain a set of four coupled equations. 
The general solution of this set of equations may be written as 



( ei(t) \ 
e 2 (t) 

V Pit) J 



6 2 



(46) 
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where Si and (3 are constants. The frequencies corresponding to small oscillations are given 
by 



W 



,2 ' 



,2 ' 



(47) 
(48) 



and 



a 2 (3 + 2w) 



<$poVZ - 6 + iu(2- 3 7 ) 



(49) 



± yj [ - 6 + a2$ ^ " + 2w - 3w 7 ] +2(3 + 2w) ( - 6 + ag$o^ " [ 3 7 - 2]) 

The static solution is stable if lo± > 0. Assuming that (3 + 2w) > we find that this 
solution is stable against anisotropic perturbations, providing that Eqs. ( ITBTOOl) are fulfilled. 
The oscillation mode which belongs to the frequency W\ is given by 



e 2 (t) 
e 3 (t) 

v m J 



/-A 



1 

V o J 



iwit 



Similarly, the oscillation mode corresponding to w 2 is: 



/ ei(f) \ 
e 2 (t) 

es(t) 



Co 



/-A 
i 

o 

V o / 



iw2t 



Finally, the oscillation modes corresponding to w± are: 



(50) 



(51) 











e 2 (t) 






e iw±t 


e 3 (t) 










V i ) 





(52) 
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where Cj (i=l,2) and C± are arbitrary constants. On the other hand, the constants A± are 
given by 



6(37 - 4) 



ag$oK' -w(2 - 3 7 ) (53) 



=F V [ - 6 + a 2 $ V " + 2w- 3wj] + 2(3 + 2w) ( - 6 + a2$ y o " [37 - 2]) 

Notice that the oscillation modes corresponding to the frequencies W\ and W2 are 
anisotropic oscillations around the equilibrium point. In these oscillations the JBD field 
remains static at its equilibrium point, $o- Cm the other hand, the oscillation modes, re- 
lated to the frequencies w± are isotropic oscillations around the same point, but where now 
the JBD field oscillate. 

We note that this stability behavior is completely different wherewith it happens with 
the ES solution, where it was found that spatially homogeneous Bianchi type IX modes 
destabilize the static solution 



19] . 



B. Scalar field 

In this case, we consider a universe dominated by a scalar field. Following a similar 
scheme to that of Sec. IHH we take a flat matter potential, Ufif), with a scalar field ^ rolling 
along its potential with a constant velocity satisfying the conditions for a static universe. We 
study the stability of this solution against anisotropic Bianchi type IX perturbation modes. 

In this case the set of Eqs. ( 13711391) becomes 

p = po + 6p(e 1 ,£ 2 , £3) ~ Po - 7o Po [si + e 2 + £3] , (54) 
P = P + <JP(ei, £ 2 , e 3 ) w P + + 7o (1 - 7o) Po) [e x + e 2 + £3] ■ (55) 



and 



7 = 7o + £2, £3) ~ 7o - 2^-^ [en + £ 2 + £3] , (56) 

Po 
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We introduce these latter expressions into Eqs. (I4"01 - (j4"3"l) and, just like above, we retain 
the linear terms in the parameters q and (3. In this way, we obtain a set of four coupled 
equations. A general solution to this set of equations could be written as Eq. 

In this case, the frequencies for small oscillation are given by 



8 



2 ' 



,2 ' 



and 



a 2 {3 + 2w) 



ajj$ V r "-2(3 + 2w) 



(57) 
(58) 

(59) 



± yK^oK,"] 2 + 4a 2 ) <W(3 + 2w) + 8w(3 + 2w) 

If io\ > the static solution is stable. Assuming that (3 + 2w) > 0, we find that this 
solution is stable against anisotropic perturbations, if the ranges expressed by expressions 
fl26l 1271) are fulfilled. Notice that these constrains are the same constrains which were found 



previously in Ref. 36J, where the stability of this static solution against homogeneous and 
isotropic perturbations was studied. 

The oscillation modes corresponding to these perturbations share similar properties than 
that the ones discussed in the previous section. In particular, they could be expressed by 
the same expressions, Eqs. (|50l IBTl 1521). but where now w±, W2 and w± are given by Eqs. ( J57l 
1581 159]) respectively, and A± is given by 



A± = — 
12 



a^o Vq +AwT\ [a§ W] 2 + 4a 2 ,$ V "(3 + 2w) + 8w(3 + 2w) 



(60) 



This modification of the stability behavior has important consequences for the emergent 
universe scenario, since it ameliorates the fine-tuning that arises from the fact that the ES 
model is an unstable saddle in GR. 

V. POLYNOMIAL JBD POTENTIAL 



As a particular but interesting example, we consider the case where the JBD potential is 
a polynomial in the scalar field 



15 



(61) 



where Co, C% and C2 are constants. Also, we consider a homogenous and isotropic closed 
universe described by a Friedmann- Robertson- Walker metric Eq. 01]). As a matter content 
we take a standard perfect fluid. 

It is interesting to notice that under these consideration and following the scheme of 



Refs. 



41 



42j we can rewrite the field equations of this model, Eqs. (JMZ]), as an autonomous 



system. In order to do so, we first rewrite Eqs. (0 [7]) together with the conservation of 
energy equation, by means of the conformal time 

dt 



a(t)' 



Thus, we obtain that 



a/_ &_y pa*_ 2w + 3 /gV C a 2 V x a 2 V 2 a 2 



(62) 



and the conservation of energy-momentum becomes 



(63) 



p' + 3- 7 p = 
a 



where, we have used the following definitions 



(64) 



Vi = Ci$. 



(65) 



and 



V 2 = C 2 $ 2 



Following Refs. 



41 



42 1 we introduce the set of variables 



(66) 
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and 



l2w + 3 $' $' 
= ^EJLIZ. = A — , (67) 
V 12 $ $ ' v ; 

a' $' 

Z» = %^ , (69) 
_ Via 2 



Now, we rewrite Eqs. (|62|) and (1631) . together with the energy- momentum conservation in 
these variables as follows 



Y + 1 



pa ~ Z Z\ Z2 
— + X 2 + — + — + — , 
3$ 3 3 3' 



(72) 



and 



X' = -2XY+[l- 3 - 1 )^+ Z ° ' Zl 



(73) 



4 7 3,4$ 3A 6A 

Differentiating Eq. (J72l) and from the equation for the X variable, together with the 
energy-momentum conservation, we obtain that 



X' = -2X Y 4 ( 1 -- 7 



F 2 + 1 - X 2 



+ 



7^0 
4 A 



7 1\ # 



+ t-7 hr+ 77-1 -T 



4 6 7 A V 4 



3A 



y = _2X 2 + f 1 - ^ 7 J [F 2 + 1 - X 2 ] + 1(Z + Z 1 + Z 2 



Z' n — 2Zri 



Z[ = 2Z 1 



X 
~A 



Y 



X 

-2A + Y 



(74) 
(75) 
(76) 
(77) 



17 



and 

Z' 2 = 2Z 2 Y . (78) 
Requiring that p > 0, we get from Eq. (172!) that 

3 3 3 V ; 

In the set of equations (1741) - (fTS!) we look for critical points. In particular, we are interested 
in critical points related to static universes which were discussed in Sect. IIHI Thus, from 
Eqs. (1741) - (fTSl) together with expression (l67j) and (1681) . the critical points correspond to 
X = Y = 0, Z = Z , Z\ = Z\ and Z 2 = Z 2 , where 

- 3 2 Z x 



-- 2 7 2' < 80 » 
and 

Z, = § - f • (81) 

Then, we have a set of critical points which represents different static universes. They 
depend on the arbitrary value of Z\. Actually, the possibility of obtaining stable or instable 
critical points depends on the value of Z\. In the following, we will give a range for the 
parameter Z\ where the corresponding solutions are stable (see Eqs. ( EH) and (15711 ) 

In order to study the nature of these critical points we linearize the set of equations 
( 17511781 near the critical points. From the study of the eigenvalues of the system we found 
that the critical points could be centers or saddles points, depending on the values of the 
parameters of the model (7 and A) and on the value of Z\. Stable static solutions correspond 
to a center, and this imposes the following conditions for the parameters A and 7, and for 
the value of Z\. 



\<1<\, (82) 



and 



q^v _ 1 o q 

^— - < Z x < -6(1 + 2A{2 + A)) + - (1 + 2Af 7 , (84) 



or 
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FIG. 1: Plot showing the evolution of two numerical solutions for a universe dominated by dust. 



4 

- < 7 < oo, (85) 
3t - 4 

< A < — , (86) 

67 - 4 ' v ; 



and 



Q'v — 1 2 9 

^— - < Z X < -6(1 - 2A{2 - A)) + ~ (1 — 2A f 7 . (87) 

These conditions are in agrement with the general stability conditions that were found 
previously in Sect |ffl] (see Eqs. (fT8l l20]) b 

In Fig. [TJit is shown a projection of the axis X and Z\. This represents the evolution of two 
numerical solutions for a universe dominated by dust. In order to satisfy the requirements 
of stability we have taken the values A = 0.008 and Z\ = —2. Here, the critical point, which 
in this graph corresponds to the point X = and Z x — Z\ — —2, represents a center. 

In Fig. [2] it is shown a projection of the axis X and Z\ of two numerical solutions for the 
case where the universe is dominated by a scalar field moving in a null scalar potential. In 
order to satisfy the requirements of stability we take A = 0.008 and Z\ = 2. As we expect, 
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FIG. 2: Plot showing the evolution of two numerical solutions for a universe dominated by a scalar 



the critical point, which in this graph correspond to the point X = and Zy = Z x = 2, it is 
a center. 

VI. CONCLUSIONS 

In this paper, we have studied the stability of the JBD static universe model. This is 
motivated by the possibility that the universe might have started out in an asymptotically 
JBD static state, in the context of the so called emergent universe models. 

We extent our previous results on stability of JBD static universe by considering spatially 
homogeneous Bianchi type IX anisotropic perturbation modes and by including more general 
perfect fluid. Contrary to GR we have found that the JBD static universe dominated by 
a standard perfect fluid could be stable against isotropic and anisotropic perturbations for 
some sort of perfect fluids, for example for dust or scalar field (inflaton). This modification of 
the stability behavior has important consequences for the emergent universe scenario, since it 
ameliorates the fine-tuning that arises from the fact that the ES model is an unstable saddle 



field. 
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in GR and prevent that small fluctuations, such as quantum fluctuations, will inevitably 
arise, forcing the universe away from its static state, thereby aborting the emergent universe 
scenario. 

In particular we found that for a standard perfect fluid with a polytropic state equation 
satisfying 7 > 2/3 it is possible to find a static solution which is stable against isotropic and 
anisotropic perturbations, with the only exception of radiation (7 = 4/3). This implies that 
for a universe dominated by dust (7 = 1), for example, we could find a solution where the 
universe is static and stable. The instability of the static universe dominated by radiation, 
although disturbing, seems not to be a problem, since in a pre-inflationary cosmological 
model, it might be possible that radiation be an element which is not dominant at all. 



361 ]. which correspond to 



Also, we found that the static JBD universe described in Ref. 
a universe dominated by a scalar field moving in a flat potential, is stable against isotropic 
and anisotropic perturbations when the JBD potential and the JBD parameter satisfy a set 
of general conditions discussed in Sect. ITVl 

Finally, we focus on a particular example of a JBD potential, which is a polynomial in 
the JBD field. This kind of JBD potential allow us to use a dynamical system approach 
for studying the stability of the JBD static universe. In this respect, we have found that 
the JBD static universe solutions are center equilibrium points. We obtained numerical 
solutions for a universe dominated by standard perfect fluids and dominated by a scalar 
field. We have considered the cases where the universe starts from an initial state close to 
the equilibrium point. The numerical solutions showed a behavior just like the expected if 
the equilibrium points are centers. 

We should stress that in this work we have studied the stability of the Jordan-Brans-Dicke 
static universe against spatially homogeneous isotropic and anisotropic perturbations, see 



Refs. 
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191 ]. Of course, also it is possible to study the stability of the JBD static universe 



against spatially inhomogeneous perturbations (scalar, vector and tensor perturbations). 
The situation for the ES solution, becomes neutrally stable against inhomogeneous scalar 
perturbations (with high enough sound speed), vector and tensor isotropic perturbations 



is 



191 ] . We expect that in the JBD case these inhomogeneous perturbations do not lead to 



additional instabilities in the same way that happens with the ES case. We intend to return 
;o this poin t in the near future by working an approach analogous to that followed in Refs. 



43, 



44 



45 



46]. 
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